Interference effects in above-threshold ionization from diatomic molecules: 
determining the internuclear separation 
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We calculate angle-resolved above-threshold ionization spectra for diatomic molecules in linearly 
polarized laser fields, employing the strong-field approximation. The interference structure resulting 
from the individual contributions of the different scattering scenarios is discussed in detail, with 
respect to the dependence on the internuclear distance and molecular orientation. We show that, 
in general, the contributions from the processes in which the electron is freed at one center and 
rescatters off the other obscure the interference maxima and minima obtained from single-center 
processes. However, around the boundary of the energy regions for which rescattering has a classical 
counterpart, such processes play a negligible role and very clear interference patterns are observed. 
In such energy regions, one is able to infer the internuclear distance from the energy difference 
between adjacent interference minima. 



I. INTRODUCTION 

The interaction of matter with an intense laser field 
(/ > 10 13 W/cm 2 ) leads to several phenomena, such as 
above-threshold ionization (ATI) or high-order harmonic 
generation (HHG). Such phenomena owe their existence 
to physical mechanisms, in which an electron reaches 
the continuum, by tunneling or multiphoton ionization, 
at an instant t' . Subsequently, it is accelerated by the 
field and driven back towards its parent ion, or molecule, 
with which it rescatters or recombines at a later time 
t Such laser-induced recombination or rescattering 
processes take place within a fraction of a laser-field cycle. 
The period of a typical near-infrared Ti:sapphire laser 
pulse is T = 2n/uj ~ 2.6fs. Thus, HHG and ATI occur on 
a time scale of hundreds of attoseconds Q ■ Hence, above- 
threshold ionization and high-order harmonic generation 
may be employed for probing, or even controlling, dy- 
namic processes with attosecond and sub-angstrom reso- 
lution. 

This fact, together with new alignment techniques, has 
opened a whole new range of possibilities for studying 
molecules in strong laser fields, employing high-energy 
photoclectrons or high-order harmonic radiation. Con- 
crete examples are the attosecond reconstruction of the 
nuclear motion in a molecule Q , the real-time imaging of 
vibrational wavepackets [4 L the tomographic reconstruc- 
tion of molecular orbitals 5| , the time-resolved measure- 
ment of intramolecular quantum-interference effects Q, 
or the determination of internuclear distances 0] • 

These applications are a direct consequence of the fact 
that a molecule possesses a very specific configuration of 
ions from which the electron may leave, or off which it 
may rescatter causing above-threshold ionization, or re- 
combine generating high-harmonics. This leads to char- 
acteristic quantum-interference patterns in the HHG or 
ATI spectra, in which structural information about the 



molecule is hidden. This is true both for polyatomic 3 
and diatomic molecules d, H, [p, ITol- Jill fl3l . 111 . [ill [la . 
El, EI El, EH, E3, El, iH Hi EI El. In particular 
for diatomic molecules, it has been shown that the high- 
order harmonic or ATI spectra exhibit overall maxima 
and minima, which are highly dependent on the spatial 
separation between both centers in the molecules, and 
can be described as the interference between two radi- 
ating point sources. In this sense, HHG or ATI by a 
diatomic molecule may be viewed as the microscopic ana- 
log of a double-slit experiment [1, El, EH ■ Furthermore, 
such features depend on the symmetry of the highest oc- 
cupied molecular orbitals, and on the alignment angle of 
the molecule with respect to the laser-field polarization 

@, 0, El, EH, El, El, El El, El, EH, El, El, E3, El, El, E3 • 

Specifically in the diatomic case, several aspects of 
this interference have been extensively studied in the 
past few years, such as the influence of the orbital sym- 
metry, the internuclear distance, the alig nment angles 
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molecular vibration 
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2a , |26| , as well as the role of the 



laser-field shape lH" 17| or polarization [l8(. Further- 
more, an adequate modeling of bound molecular states, 
in comparison with existing ionization experiments [3], 
has also raised considerable debate [H El, El, EH] ■ 

For that purpose, both the purely numerical solution 
of the time-dependent Schrodinger equation [7J, I 111. 1 121 . 
and the st rong -field approximation [ll, El, [la. Il9l.l20l. 
El, El, El, ES El, E3 have been employed. The latter 
method allows a transparent physical interpretation of 
the phenomena in question as laser-induced rescattering 
or recombination processes, and permits a clear space- 
time picture, which can be related to the classical orbits 
of an electron in a strong laser field [28|. For a diatomic 
molecule, there exist two main rescattering or recombi- 
nation scenarios: the electron born through ionization at 
the center C; upon its return may recollide and interact 
with cither the same ion (Ci), or with the other one ( Cj 
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7^ j)) Such processes have been taken into account for 
high-order harmonic generation employing a two-center 
zero-r ang e potential [lj,[2l[, using Bessel function expan- 
sions |l9| . and by means of saddle-point methods [2lll22l] . 

In this paper, we calculate the energy spectra and an- 
gular distributions of ATI produced electrons in linearly 
polarized laser fields, within the framework of the strong- 
field approximation (SFA) and the single-active electron 
approximation (SAE). We employ a zero-range poten- 
tial model similar to that in and consider both the 
direct electrons, which reach the detector without inter- 
acting with their parent molecule, and the electrons that 
suffer a single act of rescattering before reaching the de- 
tector. In the latter case, we put particular emphasis on 
interference effects: A final state with given momentum 
outside the laser field can be reached via two different 
scenarios. An electron can be born at and rescatter off 
the same center, or it can be born at one center and 
rescatter off the other. We show that the processes in- 
volving two centers, in general, obscure the interference 
patterns in the ATI spectra, in almost all energy-angle 
regions. An exception, however, is the boundary of the 
region that after tunneling is classically accessible to the 
ionized electron, in other words, the region before the 
classical cutoff. Near this boundary, the two-center pro- 
cesses yield negligible contributions, and one may identify 
very clear interference patterns. This makes it possible to 
provide a recipe to determine the internuclear distance R 
out of the angle-resolved ATI spectra. Throughout the 
article we will use the velocity gauge and atomic units 
(e = m = h = 1, c = 137). 

The paper is organized as follows: In Sec. [IT] we pro- 
vide the ATI transition amplitudes for the direct and for 
the rescattered electrons, which in Sec. IIIII arc employed 
to compute the ATI spectra. The interference patterns in 
the spectra are analyzed with respect to molecular orien- 
tation, internuclear distance, and the position of the de- 
tector with respect to the polarization of the laser field 
and the molecular axis (Sec IIII A[) . In Sec. IIIIB1 we 
present angle-resolved spectra, from which we infer the 
internuclear distance. Finally, in Sec. IIV| we summarize 
the paper. 



II. TRANSITION AMPLITUDES 

The transition amplitude for direct ionization, within 
the strong- field approximation (SFA) [29[, is given by the 
Keldysh-Faisal-Reiss amplitude [3(| 



M p = -i 



dt(*<P(t)\V\# (t)), 



(1) 



where |*o(*)) = |*o> exp(z|£7 |*)- The amplitude de- 
scribes an electron, initially in the ground state j^o}) that 
is injected in the continuum by the laser field overcom- 
ing the ionization potential \Eq\, and reaches the detector 
with final momentum p. The form of the transition am- 
plitude given here, which contains the binding potential 



V(r) rather than the interaction with the laser field, was 
first presented in Ref. [3l|. In the SFA, the final state 
with momentum p is described by a Volkov state, which 
in velocity gauge has the form 



|*W (t)> 



1 



(2tt) 3 / 2 



■ fl^ dr[p+A(r)] 2 



(2) 



In the amplitude JT]), the electron once ionized does not 
interact with the ion (that is, with the binding potential 
V) anymore. If we allow for at most one single act of 
rescattering, the amplitude (p} is replaced by 



M (o,D = _ 



dt 



d£ 



x(^(t)\VU^(t,t')V\^ (t% (3) 

Here, U^ v \t,t') denotes the Volkov time-evolution oper- 
ator, which describes the evolution of the electron in the 
presence of the external laser field, ignoring the binding 
potential. Equation §3§ incorporates direct ionization, as 
described by Eq. (JT|), as well as ionization followed by 
rescattering (for details, see, e.g., Ref. [IH). 

In order to apply Eqs. JT]) and to a diatomic 
molecule, we consider the two-center binding potential 



V(t) = Vo(r-Ri) + y (r-R 2 ), 



(4) 



where Rj (i = 1, 2) denote the coordinates of the centers 
C, (i — 1,2). For the ground-state wave function, we 
employ a linear combination of atomic orbitals (LCAO): 



<f(r) = ci* (r - Ri) + c 2 * (r - R 2 )- 
Specifically, we will use the zero-range potential 



V(v) 



—S(r)—r, 
k or 



(5) 



(6) 



whose single bound state is described by the wave func- 
tion 

!/2 1 „„ 



*o(r) 



(— 

V27T 



(7) 



with n = y/2\Eo\. The regularization operator d/dr r 
acts on the wave function to its right in order to satisfy 
the proper boundary conditions at the origin [3~H |. For 
direct ionization by a monochromatic linearly polarized 
laser field 



A(i) = Aq cos e, 



(8) 



the evaluation of the amplitude (p} is straightforward. 
Taking Ri = R/2 and R 2 = -R/2, so that R is the 
internuclear distance of the two centers, one obtains by 
expanding the exponent in the Volkov wave function in 
Eq. ([2]| into Bessel functions 



Nlj 




(9) 
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FIG. 1: (Color online) ATI spectra of the direct electrons in 
the laser polarization direction in the atomic and molecular 
case, for a molecule aligned parallel to a laser field of frequency 
u> = 0.058 a.u. and ponderomotive potential U p = 2.08 a.u. 
We consider a symmetric combination of atomic orbitals [ci = 
C2 = l/y/2 in Eq. (5)], with ionization potential |£Jo| = 0.9 
a.u. In order to facilitate the comparison, the same ionization 
potential Vb(r) was chosen in the atomic and molecular cases. 
The arrows mark the various destructive interference energies 
(n = 0,1,2) for R=6 a.u.. 



where U p = Aq/4 denotes the ponderomotive energy of 



the laser field 
to 



The prefactor, which is proportional 



is of no relevance, since we do not attempt to calculate 
total ionization rates. It is, however, worth mentioning 
that the limit of R — > is not straightforward. Below we 
will not face this limit. For a more detailed discussion, 
see, e.g., Refs. (l3l. |35|. 

The only difference between the matrix element (j9|) for 
a molecule and the corresponding matrix element for an 
atom, besides the independent prefactor, is the presence 
of the term cos(p • R/2). This term describes the inter- 
ference of electron orbits with momentum p originating 
from one or the other center of the two-center potential 
((4]). The cosine term comes from assuming a symmetric 
combination of orbitals in the ground-state wave func- 
tion ([5]) , so that ci = c% = 1/ \/2. For an antisymmetric 
combination, so that c\ = — C2, the cosine is replaced by 
a sine, leading to suppression of electrons with low mo- 
menta due to destructive interference 

HHH El- The 

interference factor cos(p ■ R/2) = cos(pi?cos6'/2) yields 
destructive interference for electrons with energies 



2 2 



(2n + 1)tt 
R cos 6 



(11) 



F = -k + exp(-fcR) /R, 



(10) 



for integer n. An illustration of the interference effect 
is given in Fig. [TJ which shows the spectrum of the di- 
rect electrons for an atom and for a symmetric diatomic 
molecule aligned parallel to the laser-field polarization. 
The clearly visible sharp dips in the spectrum, due to 
destructive interference, are indicated by the arrows in 
the figure. Next, we turn to the evaluation of the ma- 
trix element ©, which allows rescattering. With the 
two-center potential (j4|) and the symmetric ground-state 
wave function ([5]), the matrix element reads 



• - L* L* I ** I " 

x U {v \rt; r't') [V (r' + R/2) + V (r' - R/ 



For the zero-range potential ([6]), the integrations over 
space can be carried out easily (32j, which leaves a 
two-dimensional integral over the ionization time t' and 
the rescattering time t. For finite-range potentials, one 
may proceed by introducing form factors and employing 
saddle-point methods. In this case, the single-center pre- 
factors cause an overall decrease of the yield for increas- 
ing photoelectron energy. There are, however, no signif- 
icant changes in the interference patterns in comparison 
with the zero-range case, since these prefactors do not 
influence the action or the cosine factor. For a detailed 
discussion of the single-atom case, see, e.g., Ref. [36j . 

We split the eight integrals into those where the elec- 



■r e 5 rfr(p+A(r)) 2 ^ + R/2) + _ R ^ 

!)] [*o(r' + R/2) + * (r' - R/2)] e l|£o|t ' . (12) 
I 

tron rescatters off the same center from which it was ion- 
ized (r = r' = ±R/2) and those where it rescatters off 
the opposite center (r = — r' = ±R/2). We refer to the 
respective terms by M 1 ^ where i = +,— denote the cen- 
ter of ionization and j = +, — the center of rescattering. 
The integrals M ++ and Af~~, which specify the elec- 
trons coming from and rescattering off the same center, 
are essentially identical to the corresponding results for 
an atom [32|]. The structure of the molecule is reflected 
in the integrals M" 1 and M ^ which characterize the 
electrons that experience the presence of both centers. 
Evaluating the remaining two integrals over t and t', we 
substitute t' = t — r. The doubly infinite integral over t 



4 



then yields a (5-function expressing energy conservation, in terms of Bcsscl functions, wc obtain for the transition 
while the semi-infinite integration over r has to be cal- amplitudes 
culated numerically. Expanding all oscillating exponents 



M++ + M"" = 2FcoJ^)Y, S (y + u p + \ e o\~Nuj) f] J- ( 2 l+N) C^- 1 ^ 

\ ' N ^ ' l=-oo ^ 

3/2 

dr ' 



2Fcos(?-^) M^ tom \ (13) 



M + ~ 

P 



(2 \ /-oo / - \ 3/2 

x / e i^ e -i[|B |r-Ha+(2I+JV)/9-] e -^,r[l-(fiS^2)»] Jj J _ ( 



(14) 



M~ 



3/2 



(2Z+AT) 



-B4 



(15) 



The real quantities A, B± and the phases a and (3± are 
defined by 

Ae~ la = e~ 2 ^ T + — S m 2 — e~ iur , (16) 

WT 2 

R • e 

B±e~ ip± = pe± [isinuT- (1 - coswr)](17) 

T 

Upon R — > — R, we have £?± exp(— »/3±) — > 
£? T cxp(— i/?zp). Consequently, the matrix element M p 
does not change when R — ► — R. The complete matrix 
element is the sum of the terms (|T3l) - (|T5|) , 

M p = M++ + M~- + M+- + M-+. (18) 

The first two terms describe electrons originating from 
and rescattering off the same center. They are pro- 
portional to the atomic ionization amplitude M p atom ' 
[32| multiplied by the wave-function overlap F and the 
two-center interference factor cos(p • R/2), which we ob- 
served for the direct electrons in Eq. ((9]). The behavior of 
the exchange terms M+ _ and M~ + is more complicated 
and will be discussed below. 

The transition amplitude M p simplifies enormously 
when the molecule is aligned perpendicularly to the field 
so that R • e = 0. Equation (|17p shows that in this case 
B + = B- and = 0. Hence, the integrals on 



the right-hand side of Eqs. (|14j) and (jTSJ) are equal and 
Mp*" - + Mp~ + becomes proportional to 2 cos(p • R/2) just 
like M+" 1 " + M~~ . If, in addition, the electron is emitted 
perpendicularly to the field so that also p • e = 0, then 
B + = B- = and we have M p = unless N is even. 

Substituting r ^ r/w in Eqs. (fT3"|) - (fT7)) one can see 
that the amplitudes M p 3 and their sum M p depend on 
the parameters of the problem through the dimension- 
less quantities p 2 /w, U v /lu, and ujR 2 when the relative 
orientations of the vectors R, p, and e are kept fixed. 



III. PHOTOELECTRON SPECTRA 

In this section we discuss the ATI spectra computed 
employing the transition matrix elements ((T|) and (|3|), 
and a symmetric combination of equivalent centers [c\ = 
Co = 1/V2 in Eq. J5J)]. For the sake of simplicity, in 
the comparison between the atomic and molecular case 
the same ionization potential Vq is chosen. Specifically, 
we take \E \ = 0.9 a.u. in Eqs. (5) and (7) through- 
out ■ In Sec. IIII Al we perform a detailed analysis of 
the interference patterns with respect to the molecular 
orientation, rescattering scenarios, and the direction of 
electron emission, while in Sec. IIII SI we provide a recipe 
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for measuring the internuclear distance from an analysis 
of the interference patterns in the angle-resolved spectra. 



A. Analysis of the interference patterns 

As a first step, we investigate how the interference pat- 
terns are influenced by the orientation of the molecule 
with respect to the laser-polarization direction. Such re- 
sults are displayed in Figs. [2] and [3] for parallel and per- 
pendicular orientations, respectively. In both cases, we 
compare the entire ATI spectrum consisting of the di- 
rect and the rescattered electrons in the atomic and the 
molecular case. Unless stated otherwise (cf. Fig. we 
consider electron emission in the laser-polarization direc- 
tion. 

As expected from Eq. ([9|), for energies smaller than 
2U P , the main contributions to the yield come from the 
direct-ionization matrix clement ([9]). Apart from the 
interference-related factor of cos(p ■ R/2), the transition 
matrix element is identical to that obtained for a sin- 
gle atom (cf. Fig. 1). This factor is responsible for the 
sharp interference dips at the positions given by Eq. (fTTj) . 
In the plateau energy region, however, the spectra de- 
pend on the laser-field polarization in a more complex 
way, as will be discussed next. In case the molecule is 
aligned parallel to the laser-field polarization (Fig. [2]), 
the plateau is strongly enhanced in the molecular case, 
and the structure of the spectrum is very different from 
the atomic case and dependent on the internuclear dis- 
tance R. Indeed, inspection of the exchange integrals (fl4|) 
and () 1 5(1 does not reveal any simple dependence on the 
internuclear distance. Generally, for the molecular case 
there are more pathways into a given final state. For our 
case of a two-center potential, there are four pathways in 
place of one for the atomic case. If they add coherently, 
a significant enhancement can result, ideally by a factor 
of 16, which is roughly what is observed in Fig. [^be- 
fore the cutoff. The structure caused by the cosine factor 
is suppressed in the plateau region. This is caused by 
the contribution of the processes in which the electron 
is ejected from one center and rescatters off the other. 
Such processes correspond to the transition amplitudes 
Mp + and which do not exhibit the proportionality 

to the cosine that is characteristic of the the one-center 
scattering amplitudes M+ + and M p . A further par- 
ticular feature observed in this case is the displacement 
of the cut-off to higher energies with increasing internu- 
clear distance. This can be understood by the fact that 
an electron that moves from one center to the other may 
gain more energy from the field since it may be acceler- 
ated over a longer distance before it recollides. 

A strikingly different behavior is observed if the 
molecule is aligned perpendicular to the direction of the 
laser field so that R e = 0. In Fig. [3] we consider the case 
where the electron is emitted in the direction of the laser 
polarization so that p • R = 0, too. In this case, there is 
a general enhancement of the ATI yield in comparison to 



that of a single atom by roughly a factor of two for the di- 
rect electrons and a much larger factor for the rescattered 
electrons. Notice that the molecular spectrum is practi- 
cally independent of the internuclear distance [2l| . since 
the R dependence of the prefactor (fit)]) is weak for R > 2 
and the exponential of R 2 /(2t) is small for the values 
of R that we consider and the values of t that give the 
dominant contributions to the integral. The entire ATI 
spectrum does not show any interference structure, since 
the contributions from the two centers add constructively 
for p • R = 0. In fact, the cosine term in the matrix el- 
ements M°, M++ and M p simply reduces to one and 
the spectrum therefore looks like that of an atom. Specif- 
ically within the plateau, by symmetry the contributions 
from the two centers always interfere constructively. This 
results in a spectrum that is largely independent on the 
internuclear separation i?, except that the plateau is en- 
hanced, compared with the atomic case, by the existence 
of four pathways. Formally, this can be understood as 
discussed above at the end of Sec. II. 

For arbitrary p ■ R, if the electron is emitted perpen- 
dicular to the laser polarization so that p ■ e = 0, then 
we see from Eqs. (fTB| and (fTT)) that B + = B_, while 
ft- = ft+ + it. The sum of the two exchange terms then 
goes like cos(p • R/2) for even N and like sin(p • R/2) 
for odd N. This holds regardless of the orientation of 
the molecule. If the laser polarization is perpendicular to 
both the electron momentum and the internuclear axis, 
then B + = B- = 0. This implies that M p is nonzero 
only for integer TV. Each other electron peak is missing. 
Next, we discuss the interference pattern in more detail 
by analyzing the individual contributions to the transi- 
tion matrix element. In Fig. [4j we separately investigate 
the individual contributions to the amplitude (|T8j) . If 
only M+ + and M p are taken, a very pronounced min- 
inum is observed near 5U P . These matrix elements cor- 
respond to the case in which the electron is ejected from 
and rescatters off the same center, so that the minimum is 
due to the term cos(p • R/2). In the full spectrum |M p | 2 , 
however, this minimum is absent because it is filled by the 
contributions from the exchange terms M~+ and M+~ . 
For a given orientation of the molecule with respect to 
the laser field and for fixed momentum p, Fig. Q] shows 
that the contribution |M~ + | 2 of the scenario in which 
the electron is freed at the center C\ and rescatters off 
at the center C2 is different from that of \M£ | 2 where 
it is released at C2 and rescatters at C\. The same has 
been observed for hig h-order harmonic generation in a 
two-center system [13| . 

The imprints of interference can still be observed if 
the electron is emitted away from the laser polarization 
direction. This will cause, however, an overall decrease 
in the photoelectron energies for both the direct and the 
rescattered electrons. Examples are presented in Fig. [5] 
This behavior is known from atomic ionization, and its 
origin is the same in the molecular case; for a discussion, 
see, e.g., Ref. [32| . 
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FIG. 2: (Color online) Comparison of the complete ATI spec- 
tra consisting of direct and rescattered electrons in the atomic 
and the molecular case for internuclear distances of R — 2 
a.u. and R = 4 a.u.. The molecule is aligned parallel to 
the laser-polarization direction, and the electrons are emitted 
in the same direction. The arrows mark the destructive in- 
terferences (n = 0,1) of the molecule for R — 2 a.u.. The 
destructive interference for n = 1 is already in the plateau 
region, where the role of exchange terms becomes important. 
The remaining parameters are the same as in Fig. 1. 
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energy in units of U 



FIG. 3: (Color online) The same as Fig. [2] but with the 
molecular axis perpendicular to the laser polarization. The 
electrons are emitted parallel to the laser polarization. 



B. Determining the internuclear distance 




energy in units of U 



FIG. 4: (Color online) Individual contributions of the var- 
ious rescattering scenarios to the total amplitude (|18|l for 
a diatomic molecule with internuclear distance R = 2 a.u. 
aligned parallel to the laser-field polarization, for the same 
molecular and field parameters as in Fig. [2] The electrons are 
emitted in the polarization direction. The arrows mark the 
respective cutoff energies for the various transition amplitude 
matrix elements. The inset at the lower left is an enlargement 
of the region near the cutoff where the direct terms and the 
exchange terms differ in a characteristic fashion, allowing for 
the determination of the internuclear separation. 




FIG. 5: (Color online) Electron yield for the parameters of 
Fig. 1 and internuclear distance R — 2 a.u. for different 
emission angles ip with respect to the polarization of the laser 
field. The molecule is aligned parallel (perpendicular) to the 
laser field in the upper (lower) panel. 





FIG. 6: (Color online) Angle-resolved ATI spectra on a log- 
arithmic scale for a diatomic molecule with internuclear dis- 
tances R = 2 a.u. and R = 3 a.u. (middle and bottom 
panels, respectively), aligned parallel to the laser-field polar- 
ization, compared to the single-atom case (upper panel). The 
binding energy is Eo = 0.9 a.u. in all cases, and the laser fre- 
quency and the ponderomotive potential are to = 0.058 a.u. 
and Up — 2.08 a.u., respectively. The plotted lines depict the 
minima (solid lines) and the maxima (dashed lines) of the 
energy distribution given by Eq. (flTj) . 



FIG. 7: (Color online) Angle-resolved ATI spectra on a log- 
arithmic scale for a diatomic molecule with ionization poten- 
tial Eo = 0.9 a.u. and internuclear distances R = 4 a.u., 
R = 6 a.u. and R = 8 a.u. (upper, middle and bottom 
panels, respectively), aligned parallel to the laser-field polar- 
ization. The field parameters are the same as in the previous 
figure. The plotted lines depict the minima (solid lines) and 
the maxima (dashed lines) of the energy distribution given by 
Eq. CLD- 



For a complete picture of the angle-resolved ATI spec- 
trum, not restricted to emission in particular directions, 
we will now present density plots. While they invariably 
imply loss of fine details and depend on the positioning 
and gradient of the false-color scale, they give a compre- 
hensive overview of the general structure. We restrict 
ourselves to the case of parallel alignment [4l[. In this 
case the spectrum is symmetrical with respect to the in- 
ternuclear axis. It is obvious that the electrons with max- 
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FIG. 8: (Color online) Enlargement of a limited energy-angle 
region of Fig. [7] for R = 8 a.u. (lowest panel) with increased 
resolution. The indents of Fig. [7] are distinctly visible as val- 
leys deeply cut into the high ridge that precedes the cutoff. 



imal kinetic energy will be detected in the direction of the 
laser field. 

The angle-resolved spectra displayed in Figs. [6] and [7] 
are very intricate and do not exhibit any simple struc- 
tures. They depend strongly on the internuclear sep- 
aration but do not, on a first inspection, lend them- 
selves in any obvious way to the assignment of a specific 
value of R to a given spectrum. Especially, owing to the 
presence and magnitude of the exchange terms (|14p and 
(fl"5|) . the two-center interference, which is expressed in 
the cos(p • R/2) term, is not immediately visible. How- 
ever, looking more closely, one can observe a very distinct 
manifestation of this term just near the classical bound- 
ary of the spectrum. Roughly, the latter agrees with the 
boundaries of the colored areas in the various panels of 
Figs. [5] and [7J We observe well-defined indents in the 
overall smooth curve that defines the classical boundary. 
The positions of these indents and, especially, their sepa- 
rations agree quite well with the interference minima pre- 
dicted by Eq. (fTTj) . The figures show that the separation 
5E between the indents (on the scale of U p ) monotoni- 
cally decreases with increasing R. Hence, by comparing a 
measured angle- resolved spectrum with Figs. [5] and [7J we 
can infer the internuclear separation. For the parameters 
underlying Figs. [5] and \7\ the resulting function R(5E) is 
given in Table fl] Fig. [8] exhibits an enlargement of the 
relevant area around the classical cutoff for the case of 
R = 8 a.u. with a higher resolution of the electron yield. 
The indents are very clearly visible like valleys that cut 
into the drop of a plateau on a topographical map. 

An analytical formula for R(5E) can be gained from 
an analytical formula for the classical cutoff energy E(6) 
as a function of the angle 9. Intersecting this with the 
energies of the interference minima given by Eq. (jTTJ) al- 
lows one to determine the function R(SE) in dependence 
of the parameters of the problem. For the case of an 
atom, such a formula for E{0) is actually known [3{|. At 
least for R < 6 a.u., Figs. [S] and [7J show that this classical 
boundary does not depend very strongly on the inter- 



nuclear separation R, so that the atomic result could be 
employed. However, even with this simplification, the re- 
sulting formula is quite complicated and we refrain from 
presenting it here. 

The question arises of why near the classical bound- 
ary the interference term cos(p • R/2) roughly multiplies 
the angle-resolved spectrum, like it does for the direct 
electrons. The answer can be inferred from Fig. 2] The 
total ionization amplitude M p is the superposition (fTS"]) 
of four different scenarios such that the electron starts 
from and rescatters off one or the other center. The two 
contributions (fT3"|) where they start from and rescatter 
off the same center are identical except for the geomet- 
rical phase, which leads to the cosine in Eq. (fTS")) . In 
contrast, the other two contributions (fF4")l and (fTS)) are 
uncorrelatcd since they are generated by geometrically 
different scenarios. Their magnitudes are different and 
almost nowhere do they exhibit a significant construc- 
tive interference. The two contributions p3[) arc individ- 
ually large when the long orbit and the short orbit add 
constructively, as is the case specifically just before the 
classical cutoff. In this case, they dominate the other 
two terms (fF4"f and (fT3|) by a factor of the order of 2 to 
4. Hence the complete spectrum distinctly exhibits the 
geometrical interference, which is expressed in the factor 
cos(p • R/2). 



Internuclear distances 
of the molecule 


Energy differences of the indents 
at the spectral boundary 


R = 2 a.u. 


SE& 


4.5U P 


R = 3 a.u. 


SE^ 


3.2U P 


R = 4 a.u. 


8E^ 


2AU P 


R = 6 a.u. 


SE^ 


1.3U P 


R = 8 a.u. 


5E^ 


1.1U P 



TABLE I: Energy differences between adjacent indents 
around the classical boundary of the angle-resolved spectrum 
of Fig. [6] The differences are taken, for each internuclear 
separation, by starting with the first indent for f > 0° as a 
function of the energy. 



IV. CONCLUSIONS 

We have analyzed ATI spectra for a two-center 
molecule in a linearly polarized laser field. The terms 
of the two-center wave function contributing to the in- 
terference structure within the SFA formalism could be 
identified as well as the absence of the interference struc- 
ture throughout most of the plateau region. We have 
shown that the angle-resolved spectra can be used to de- 
termine the internuclear distance of a molecule aligned 
with the laser field, by reading off the energy differences 
between subsequent interference minima at the classical 
boundary of the spectrum. 

The validity of this method depends upon how close 
to reality is the angle-resolved spectrum calculated for 
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our model molecule. Certainly, the spectrum of the di- 
rect electrons cannot be trusted for this purpose. How- 
ever, high-order ATI of an atom is well described by the 
SFA and a zero-range potential, especially near the classi- 
cal boundary [32j ; for a comparison of spectra calculated 
from the SFA with the solution of the time-dependent 
Schrodinger equation, see Ref. [40j for the case of an 
atom. Experimentally, application of the method re- 
quires a high detection efficiency that allows one to ob- 
tain a sufficient number of counts down to the classical 
boundary. 

The problem of how to extract the internuclear sepa- 
ration from a diffraction pattern has been addressed by 
a different method in Ref. Q, employing the numerical 
solution of the time-dependent Schrdinger equation. In 
Q, however, it is necessary to compute a radial distribu- 
tion function from the diffraction intensity, whereas, with 
the method discussed in this paper, one may determine 
the internuclear distances directly from the spectra. The 
method suggested in Ref. Q has the advantage that it 
analyzes direct electrons and, therefore, does not require 
exceptionally high detection efficiency. 

Finally, in a real physical system, there exist additional 
effects, which have not been incorporated in this model 
and may alter the interference patterns. Molecular vi- 
bration, for instance, causes an intensity loss in the high- 
harmonic signal [25| , which may lead to a blurring in the 
patterns. However, recently, numerical ATI computa- 
tions in which such an effect is included have shown that 



for Hj the angle-dependent interference patterns related 
to the double-slit physical picture remain distinguishable 
in the case considered [131 • Generally, the amount of 
blurring depends on the rigidity of the vibrational poten- 
tial. Since the period of vibrations is much longer than 
the laser period, with a few-cycle laser pulse our method 
could be used to track a vibrational wave packet or the 
dissociation of a molecule Q . Another feature which has 
not been incorporated in our model is the dependence of 
the ionization potential on the internuclear distance. In 
fact, we have taken Eg to be constant, whereas, in reality, 
it decreases with R This feature, however, will only 
cause an overall energy shift in the spectra. Therefore, 
it will not affect the distance between two consecutive 
minima or maxima in the patterns for constant internu- 
clear distance (Figs. 6 and 7). Therefore, we expect our 
method to be applicable to real physical systems and a 
wide parameter range. 
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